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ABSTRACT 


The capability of a submersible vehicle to accurately maintain its commanded depth, 
in a variety of operating speeds, depths, and loading conditions is critical for mission 
accomplishment. Below a certain critical speed, a phenomenon known as dive plane 
reversal occurs, where depth response changes sign with respect to a given dive plane 
command. This thesis builds on previous studies of the phenomenon and it presents a 
comprehensive study of dive plane reversal envelopes in the presence of external forces 
and moments on the vehicle. Based on these results, rational design and operational 


decisions can be made in order to avoid unpredictable vehicle response. 
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I. INTRODUCTION 


In submarine operations, one attribute critical to mission accomplishment is a 
submarine’s ability to accurately maintain commanded depth. This is especially 
important when operating in shallow waters or at periscope depth when failure to 
maintain ordered depth could lead to easy detection of the sub or grounding. 

At low speeds, submarines display a phenomenon commonly referred to as dive plane 
reversal. When operating at speeds sufficiently slow, dive plane commands must be 
reversed in order to achieve the desired behavior in depth change. The speed where this 
transition occurs is known as the critical speed. 

Control systems for submersible vehicles must have the ability to recognize when 
the vehicle is operating at about or below the critical speed. Failure to do so could lead to 
erratic behavior of the vehicle. This is especially important when designing the control 
system for an unmanned submersible, which operates with little or no human 
intervention. In this case, recognition of the critical speed as well as dive commands once 
below this threshold must all take place automatically. 

Using previous studies on the phenomenon of dive plane reversal [Refs 1&2], this 
thesis presents a comprehensive sensitivity study of dive plane reversal envelopes under 
the influence of a variety of external forces and moments on the vehicle. Submersible 
control systems utilizing the results displayed in this thesis will be better equipped to 


maintain ordered depth at low speeds and therefore avoid unpredictable vehicle response. 








Il. VEHICLE MODELING 


A. INTRODUCTION 


For the purposes of this study in relation to the maneuvering and kinematic 
analysis of the vehicle, the following assumptions apply: 

1) The vehicle behaves as a rigid body. 

2) In regards to the acceleration of the center of mass of the vehicle, the Earth’s 
rotation is negligible. _ 

3) The primary forces that act on the vehicle are inertia, gravity, hydrostatics and 


hydrodynamics. 
B. COORDINATE SYSTEMS AND POSITIONAL DEFINITIONS 


In order to study the motion of a submersible vehicle a system for describing 
its position must first be defined. A global coordinate frame, OXYZ, is defined with 
origin O, and a set of axes aligned with directions North, East and Down. Aligned with 
this is a right-hand reference frame with unit vectors J, J and K. The necessity that 
the K direction is chosen as positive downwards results from the convention of 
referencing a submarine’s depth downwards from the water’s surface. Neglecting the 
earth’s rotation, the IJK coordinate frame becomes an inertial reference frame in which 


Newton’s Laws of Motion will be valid. A vehicle’s position, R,, in this frame will have 
the vector components, R, =[X,1+Y,J +Z,K]. From the perspective of the vehicle, a 


standard convention will be used that places the Y axis to the right while facing along the 
X axis, and the Z axis is positive downwards. 

A body centered coordinate frame O’,,,, is defined with the origin O’, located on 
the centerline of the vehicle, translating and rotating with the vehicle. The vehicle’s 


center of mass may or may not be concurrent with the origin O’. This origin will be the 





point about which all vehicle body forces will be computed. By convention, the body 


centered coordinate frame will have its origin O’ at the ship’s center, with unit vectors, i, 


j and k. The x axis is defined parallel to the ship’s longitudinal axis with the z axis 
vertically down. Generally, the ship’s center of gravity and center of buoyancy are not 
located at the body centered origin, nor are they collocated. These points are defined by 
Cg and Cg respectively. The locations of the center of buoyancy and center of mass have 
considerable bearing on this study. In Newtonian physics, the rate of change of linear 
momentum of the center of mass and moments about the body center of mass is equal to 
the rate of change of angular momentum. This is significant as the center of buoyancy of 
a submerged body is defined by its shape, while its center of gravity is determined by the 


distribution of weight over the body. 





Figure 2-1 Coordinate Axes Representation 








C. ANGULAR POSITION IN THE GLOBAL COORDINATE FRAME 


For the purposes of this thesis, the attitude of the vehicle is described by three 
angles relative to the global coordinate frame. The three terms used to describe attitude 
are azimuth, elevation and spin. Azimuth is the direction in the horizontal plane that the 
vehicle points. Elevation is the angle above or below the horizontal plane that the vehicle 
points. Finally, roll is the angle that the vehicle is rotated about its own longitudinal aXIS. 

It should be noted that when a vehicle approaches a purely vertical angle of attack 
(ie an elevation angle of 90 degrees), spin angle cannot be differentiated from azimuth 
angle. For computer simulations and control systems, a more complicated system to 
describe angular position is usually required. However, for the purposes of this study, 


azimuth, elevation and spin are sufficient to describe a vehicle’s attitude. 
D. ROTATIONAL TRANSFORMATIONS 


When describing positional, velocity and acceleration vectors, it is frequently 
more useful to describe them in relation to a body rather than that of the global coordinate 
frame. To facilitate this, a transformation system must be defined to relate global 
coordinates to body centered coordinates. The forward transformation starts with the 
global coordinate frame. Three rotations are performed to align a frame parallel to the 
body centered coordinate frame. Azimuth rotation, wy, is defined as a positive rotation 
about the global Z axis. Elevation, 0, is defined as positive rotation about the new Y axis. 
Finally, roll, o, is defined as positive rotation spout the new X axis. These three rotations 
uniquely define the angular orientation of the body. 


Consider a position vector, R, , defined in the global coordinate frame by 
R, = [xo, yo, Zo ]’. It can be seen that R, will have different coordinates when related to 


a body centered coordinate system that has been rotated by an azimuth angle of y about 


the global Z axis. 





Referring to Figure 2.2, the position vector, R,, is related to the global 
coordinate frame using the following relationships: 
Xj} = Xp cOSW + yo siny (2.1) 


yr = -Xo SINW+ yo CosW (2.2) 


Xo 





Figure 2-2 Azimuth Rotation 


These relationships can be represented in matrix form using the rotation matrix operation, 

RK, = [1,2] Ro; (2.3) 
where [T,,, ] is an orthogonal transformation{Ref 1]. Premultiplying any vector, Res by 
the transformation [7,,, ], yields the components of the vector in the rotated coordinate 


frame. Extending this rotation to three directions yields a combined total rotation 


transformation, 


T(8,0,w) = T(®)T(O)T(y). (2.4) 


Expanding equation 2.4 results in the following matrix: 











cosy cos@ sinycos@ —sin@ 
cosysin@sing—-sinycos@ sinwsin@sing cos@sing 
cosysin@cos@+sinysing sinwsin@cosd cos@cos@ 


Using the above transformation, any vector in the global coordinate frame can be 
expressed in coordinates of a rotated reference frame using 


Rtas) [T(8,0,y)] Regs : (2.5) 


E. KINEMATICS 


In order to study the motion of an object, definitions of translational and 
rotational velocities must be established. These definitions are utilized when further 
studying the effect of forces and accelerations on the object. Translational velocity is 


defined as, 
xX | 
R=|Y (2.6) 

Z 


Using the rotational transformations described in the previous section, this global velocity 
vector can be translated into a body fixed reference frame. Let [u, v, w]’ be the body 
fixed components of the velocity vector, where u is the vehicles forward speed, v is its 
side slip velocity (sway) and w is its speed in relation to its local vertical axis (heave). 


The body fixed velocity can be determined via the forward transformation: 


u X 
v|=T(¢,0,w)| Y (2.7) 
w Z 





Transformation back to global coordinates 1s accomplished by the reverse transformation: 


X u 
Y |=T'(¢,6,w)| v (2.8) 
Z w 


Using these relations, the global motion of a vehicle can be described in terms of its 
angular attitude and its body centered velocity components. 

Transformation of rotational velocities from global coordinates to body centered 
coordinates is more complicated. Because elevation angle, 8, is made about an 
intermediate y axis based upon azimuth and spin angle, 0, is made about an intermediate 
x axis based upon the elevation angle, the transformation is a three step process. Let @ be 
the rotation rate of the body relative to its own body aligned reference frame and p, q and 
r be the azimuth, elevation and spin components of @. @ is calculated from global 


angular velocities via the following equation: 


p 0 0 d 
q|=ZXTOVEY) 0|+TO@TE} 4 |+T@} 0 (2.9) 
r W 0 


Simplifying equation 2.9 yields: 


p —wsinO+¢ 
qi=| wsin@+@cosd |. (2.10) 
r wcosOcos¢—-Osing 


Calculating global angular velocities from body centered angular velocities is 


accomplished by inverting equation 2.10: 








d p+qsin@gtan@+rcos@tand 
= gcosé—rsing (2.11) 
(gsing@+rcos@)/cosé 


<. @. 


F. EQUATIONS OF MOTION 


Newtonian physics is the basis for describing motion of an object and relative 
forces involved. In order to apply F = mA to translational motion, acceleration of an 
object within a rotating reference frame must first be defined. This acceleration is 


derived by differentiating the position vector R, with respect to time twice. 


Ro =V+tOxp,et+t@xOxp,+@xy. _ (2.12) 


Forces applied to an object are equated to the product of this acceleration and mass of the 


object. 
F =m +®xp, t+WOxOxp, +OxD, } (2.13) 


For a submersible vehicle, this force vector is comprised of gravitational, hydrostatic and 
hydrodynamic forces acting on the vehicle. 

For rotational equations of motion the rate of change of angular momentum about 
an object’s center of mass is equated by the sum of moments applied about the center of 


mass. The mass moments of inertia for an object are defined as follows: 


by xy i 
ct a oe (2.14) 
i We, te 


Angular momentum for an object 1s defined as follows: 


H, = 1,0. (2.15) 


M, =H, +7qx(nk, | (2.16) 


Using equations 2.12 and 2.16, the total applied moments on an object about the origin is 


given as follows: 


H, =1,0+@0xH,. (2.17) 


where the acceleration of the global position vector 1s given by, 


R, =V+@xv. . (2.18) 


Combining these elements, the rotational equation of motion becomes: 


M, =1,0+@x(I,@)+m\p, xv +P, xox} (2.19) 
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G. VERTICAL FORCES 


The two vertical forces always present in submersible dynamics are submersible 
weight and buoyancy. For a submersible to be neutrally buoyant, the magnitudes of its 
weight and buoyancy must be equal. A submersible’s weight is defined by its mass and 
acts through the center of gravity of the submersible. A submersible’s buoyancy is 
defined by its shape and acts through the submersible’s center of buoyancy. Usually, but 
not always, the centers of buoyancy and gravity are located on the same vertical line of 
action. For a submersible to be statically stable, the center of buoyancy must be located 
above the center of gravity. As such there is usually a vertical separation between the 
centers of buoyancy and gravity. 

Assume that a submersible is operating at a constant depth at a pitch angle 8 = 0 
and roll angle = O and that its centers of buoyancy and gravity are act through the same 
vertical line of action. Any change in pitch or roll will cause a lateral separation of the 
respective lines of action of buoyancy and gravity. As a result, a moment will develop | 
acting to force the submersible back to its original orientation. The resolved vertical 
force based upon weight, buoyancy, pitch and roll is, 

~sin@ | 
f, =(W -B)| cos@sing (2.20) 
cos 8 cos@ 
and always acts in the global vertical direction. Note that if the magnitudes of weight and 
buoyancy are equal, this resolved force is zero. The resolved moment based upon weight, 


buoyancy, pitch and roll is, 


—sin§ ~sin@ 
m, =Wp, X| cos@sing |— Bp, X| cos@sing |. (2.21) 
cos @cos@ cos @cos@ 


Because there is a vertical separation between the centers of gravity and buoyancy, a 
moment will develop with changes in pitch and roll even if the magnitudes of weight and 


buoyancy are equal. 


1] 
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Ii. DIVE PLANE REVERSAL 


In order to study the phenomenon of dive plane reversal, dive plane function must 
first be understood. There are several ways to control the depth of a submersible vehicle. 
Buoyancy can be altered or adjustments can be made to either the fairwater or dive 


planes. The primary method of depth control is the utilization of dive planes. 
“I 
Fero 


Fe 


Figure 3-1 Submersible’s initial configuration 


Assume a submersible vehicle is travelling through the water with a speed v ata 
constant depth with a pitch angle of 6 = 0. Referring to Figure 3-1, the forces acting on 
the submersible are its buoyancy, ‘gravity and its propulsion force. To descend, a 
submersible’s dive planes are deflected to an upward angle of attack a. An upward force, 
F 


plane 


which is a function of & and v* results on the dive planes. Acting through the 


submersible’s metacenter, F nlane Creates a moment, M which acts to pitch the 


lane plane 


submersible forward to some negative pitch 0. Referring to Figure 3-2, as the submersible 
pitches forward, the centers of buoyancy and gravity which were vertically aligned 
originally achieve a small amount of longitudinal separation and a couple develops. This 


moment, M rich acts to oppose the action of M and is strictly a function of the 


plane 
submersible’s pitch angle 8. Eventually, as the pitch angle increases, the magnitudes of 
M iter ANd M___,,. Will equalize and the submersible will stabilize at that pitch angle. 


plane 
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Mpilane 


Fpiin 


Mpitch 


Figure 3-2 Submersible’s Diving Configuration 


Dive plane reversal is a phenomenon that occurs below a speed known as the 


critical speed. Below the critical speed the hydrodynamic force F., and subsequent. 


plane 


moment M nlane aren't large enough appreciably alter the pitch of the submersible. 


Consequently, the only force acting to change the depth of the submersible is F A 


plane * 
submersible attempting to dive using its dive planes below the critical speed will actually 
rise towards the surface. 

Onboard the Navy’s operational submarine fleet, dive plane reversal is well 
understood. When operating below a submarine’s critical speed, its operators know to 
take different actions to control depth. However, with unmanned submersible vehicles, 
the phenomenon of dive plane reversal leads to complicated problems in the area of 
automated control systems. An unmanned submersible’s control system must be able to 
recognize when the submersible is operating below the critical speed. Additionally, it 
must take appropriate action to control and maintain depth. Without these important 
considerations, control actions could lead to erroneous results and instabilities in vehicle 


control. For a vehicle operating in restricted waters, shallow depths and engaged in 
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delicate activities such as mine hunting or sweepings, the importance of a control system 


to be stable below the critical speed increases drastically. 


I> 
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IV. RESEARCH 


Under certain conditions, submersible vehicles can be influenced by external 
forces and moments due to operating conditions. Surface effects are vertical forces 
resulting from a pressure differential as a result of operating near the surface. 
Additionally, external moments may be produced due to surface effects depending upon 
the angle the vehicle’s longitudinal axis makes with the surface. Similarly, bottom effects 
are vertical forces resulting from a pressure differential as a result of operating near the 
ocean floor. 

The focus of this thesis is to explore the effects that external forces and moments 
have on a submersible’s critical speed. This research utilized the MATLAB program 
forcer.m. This program numerically calculates the critical speed envelopes of the 
DARPA SUBOFF model based on a set of input characteristics. The DARPA SUBOFF 
model is used in this simulation as a generic body of revolution for which a set of 
hydrodynamic coefficients and geometric properties is readily available. The program 
forcer.m is a modified version of bifsu.m [Ref 1] and was modified to utilize user defined 
inputs of external vertical forces and moments. During simulation the effect of these 
external forces and moments can be observed. | | 

For this study the magnitudes of weight and buoyancy were assumed to be equal. 
Using body centered coordinate frame aligned with the submersible’s geometric center 
Newton’s equations of motion are expressed as: 

m(w—Uq-26q° —X%¢4) =Z,q+Z,w+ZUq+Z,Uw 
aC. | 10) Eas + F, cos@ +U°Z,6 ile 


stem 


lw 


I ,qt+mz,wq—mx,(w-Uq)=M ,q+M,wt+M Uq+M,Uw+x,F, cos6 
(4.2) 


Ww-X 


bow 
—(zgW -z,B)sin9+U°M,5-C, | b(x) 


stem 


3 
) xdx —(x,W — x,B)cos@ 
lw — xq| | 
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In equations (4.1) and (4.2), mis the vehicle’s mass, U is the vehicles forward speed, g is 
the pitch rate or 6, Z is the heave force, w is the heave velocity I, is the vehicle’s mass 
moment of inertia and 6 denotes the stern plane deflection. The coordinates for the 
vehicle’s centers of gravity and buoyancy are denoted by (%G ,zg ) and (xg io ) 
respectively. b(x) denotes the beam of the vehicle as a function of length. To simplify 
analysis, coefficient of drag, Cp is assumed to be constant over the length of the vehicle 
[Ref. 1]. The external vertical force and point of application are expressed by F, and xr. 
Subscripts represent the derivative of the variable with respect to the subscripted value. 
For example, M, represents the derivative of pitch moment with respect to heave 
velocity. 


Using these conventions, the rate of change of depth is calculated by: 
z=—-Usin6+weos6. — (4.3) 


To simplify analysis, the forward velocity of the vehicle U is assumed to remain constant. 


Assuming steady state, equations (4.1), (4.2) and (4.3) algebraically reduce to: 


(Z,M;-M,,Z;)Uw-C,A,(M5—x,Zs5)ww+LF.M, 
+(xpW -x,B+F,x,)Z;]cos@+Z;(zoW —z,B)sin@ =0 


(4.4) 


ence: 
WwW U 
sin @ = ———_——-. , cos 8 = 
Vw? +U? vw? +U? 
1 
A. = | b(x)dx Xx, = re aia 


The effects of control surface saturation are not taken into account in equation 
(4.4). As a result, steady state solutions in pitch are allowed to achieve unrealistically 


high values. In order to account for contro] surface saturation, 6 = 6_,, is utilized, where 


f 








6, =1+0.4 radians. As aresult, w#U tan6. Steady state values are calculated from a 


heave velocity equation: 
Z,,Uw-C,A,ww - F,cos6+Z;U°6,,, =0 (4.5) 


and a pitch rate equation: 
M ,Uw-C,x,A,U°wlw|-(x.W -x,B-x,F,)cosé (4.6) 
—(z,W -—z,B)sin@+M ,U’6,,, =0 

which are utilized in the calculation. 

During calculation, forcer.m identifies a value of Xgb/L corresponding to the 
critical speed based on loading conditions. Xgb/L is simply a nondimensional parameter 
describing the longitudinal separation of the submersible’s center’s of gravity and 
buoyancy divided by its overall length. Critical speed is converted to a Froude number 


defined as, 

| Fn =U 1/4) 826, - (4.7) 
Calculations are performed over a range of speeds. Critical Froude number is 
subsequently defined as, 


1/2 
Fn ‘ao | (4.8) 
(M Z;—-Z M 5)g 


Although Froude number is strictly associated with surface vessels and does not normally 
apply to fully submerged bodies, it is applicable here. Froude number is a 
nondimensional term describing the ratio of hydrodynamic and hydrostatic forces. In the 
case of equation (4.4), that is precisely what it is quantifying. A submersible’s critical 
speed will depend primarily on its metacentric height. Another advantage of using 
Froude number is that because it is a nondimensional term, it can be appropriately scaled 
for use with different submersibles. After the program has performed its calculations, it 


plots critical Froude number as a function of Xgb/L. 
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To research the affect of external forces and moments on critical speed, forcer.m 
is run repeatedly. By varying force magnitude, force direction, point of application 
coefficient of drag, conclusions can be made as to the influence these parameters have on 


the determination of critical speed. 
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V. RESULTS 


The simulation was originally run with no external forces and moments present to 
produce a baseline plot of the critical speed envelopes. Figure 5-1 displays these 
envelopes at three different coefficients of drag. Referring to Figure 5-1, the peaks of the 
plots are all centered at an Xgb/L value of 0. As the longitudinal separation of the centers 
of buoyancy and mass increases the corresponding critical speed decreases. Similarly, as 
the coefficient of drag increases, the critical speed decreases. 

Figure 5-2 displays the critical speed envelopes in the presence of a upwards 
vertical force at a coefficient of drag of 0.0. These forces are applied at the submersible’s 
metacenter and hence incur no moment on the submersible. Referring to Figure 5-2, 
increasing in force magnitude yields a shift of the critical speed curve to a positive value 
of Xgb/L as well as a slight decrease in the peak value of critical speed. The relative 
magnitude of these forces to the submersible’s overall weight is very small (< 0.5%). 
Yet, even at these small magnitudes, significant shifts in the critical speed envelopes 
result. Similar trends are displayed on Figures 5-3 and 5-4 for coefficients of drag of 0.1 
and 0.3 respectively. 

Figure 5-5 displays the critical speed envelopes in the presence of a downwards 
vertical force. Similarly as before with upwards forces, there is a decrease in critical 
speed with increasing applied force. However, with negative vertical forces present, the 
shift of the critical speed curve is to a negative value of Xgb/L. Similar trends are 
displayed on Figures 5-6 and 5-7 for coefficients of drag of 0.1 and 0.3 respectively. 

Figure 5-8 displays the critical speed envelopes in the presence of a positive 
external moments at a Cp of 0.0. For comparison, these moments were produced by 
forces equal in magnitude as in previous plots but applied at a point a quarter of the 
submersible’s length along its longitudinal axis. Again, a trend of decreasing critical 
speed is observed in the presence of greater external force. However, the Xgb/L shift is 
less pronounced when comparing external moments to external force. Similar trends are 


displayed on Figures 5-9 and 5-10 for coefficients of drag of 0.1 and 0.3 respectively. 
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Figure 5-11 displays the critical speed envelopes in the presence of a positive 
external moments at a Cp of 0.0. Predictably, critical speed decreases as the magnitude of 
the moment is increased. The Xgb/L shift is negative and not as pronounced as that in the 
presence of external negative vertical force. Similar trends are displayed on Figures 5-12 


and 5-13 for coefficients of drag of 0.1 and 0.3 respectively. 
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Critical Speed Envelopes for Various Cd with no External Vertical Forces 
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Figure 5-1 Baseline Critical Speed Envelopes 
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Critical Speed Envelopes for Cd = 0.0 at Various Vertical Forces _ 
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Figure 5-2 Critical Speed Envelopes Subject to Various Positive Vertical Forces at Cp = 0.0 
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Critical Speed Envelopes for Cd = 0.1 at Various Positive Vertical Forces 
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Figure 5-3 Critical Speed Envelopes Subject to Various Positive Vertical Forces at Cp = 0.1 
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Critical Speed Envelopes at Cd = 0.3 at Various Positive Vertical Forces 
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Figure 5-4 Critical Speed Envelopes Subject to Various Positive Vertical Forces at Cp = 0.3 
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Critical Speed Envelopes for Cd = 0.0 at Various Negative Vertical Forces 





Figure 5-5 Critical Speed Envelopes Subject to 


Various Negative Vertical Forces at Cp = 0.0 
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Critical Speed Envelopes for Cd = 0.1 at Various Negative Vertical Forces 
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Figure 5-6 Critical Speed Envelopes Subject to 


Various Negative Vertical Forces at Cp = 0.1 
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Critical Speed Envelope Plot for Cd = 0.3 at Various Negative Vertical Forces 
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Figure 5-7 Critical Speed Envelopes Subject to 


Various Negative Vertical Forces at Cp = 0.3 


2g 





Critical Speeds for Cd = 0.0 for Positive Vertical Forces at L/4 forward of amidships 
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Figure 5-8 Critical Speed Envelopes Subject to 


Various Positive Moments at Cp = 0.0 
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Critical Speeds for Cd = 0.1 for Positive Vertical Forces at L/4 forward of amidships 
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Figure 5-9 Critical Speed Envelopes Subject to 


Various Positive Forces at Cp = 0.1 
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Figure 5-10 Critical Speed Envelopes Subject to 


Various Positive Moments at Cp = 0.3 
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Critical Speeds for Cd = 0.3 for Positive Vertical Forces at L/4 forward of amidships 
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Critical Speeds for Cd = 0.0 for Negative Vertical Forces at L/4 forward of amidships 
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Figure 5-11 Critical Speed Envelopes Subject to 


Various Negative Moments at Cp = 0.0 
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Critical Speeds for Cd = 0.1 for Negative Vertical Forces at L/4 forward of amidships 
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Figure 5-12 Critical Speed Envelopes Subject to 


Various Negative Moments at Cp = 0.1 
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Critical Speeds for Cd = 0.3 for Negative Vertical Forces at L/4 forward of amidships 
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Figure 5-13 Critical Speed Envelopes Subject to 


Various Negative Vertical Forces at Cp = 0.3 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The main conclusions of this work can be summarized as follows: 

1) Positive vertical forces tend to lower a submersible’s critical speed while 
shifting the critical speed envelope to the right on the Xgb/L axis. 

2) Negative vertical forces tend to lower a submersible’s critical speed while 
shifting the critical speed envelope to the left on the Xgb/L axis. 

3) Positive and negative moments tend to lower a submersible’s critical 
speed. The critical speed envelope shift is similar to that of positive and 
negative external forces. However, the shift is less pronounced. 

The most notable feature of all results that were generated in this work is the 
extreme sensitivity of the parameter envelopes to the magnitude of the external forces on 
the vehicle. The regions where dive plane reversal is predicted to occur seem to shift 
dramatically even under the smallest external excitation. Although such a shift can be — 
potentially very difficult to deal with in an automatic control system, it can be utilized to 
the designer’s advantage since, in principle, it could establish innovative ways of accurate 
determination of external forces and moments. Such an identification scheme could prove 
very useful in shallow water mine hunting and other similar operations in littoral waters 
or the surf zone. Such advanced disturbance observation schemes along with robust 


control studies are recommendations for further research. 


B. RECOMMENDATIONS 

Utilize the findings of this thesis in the development of robust control systems for 
unmanned submersible vehicles that are operating at low speeds in shallow waters. Such 
control systems must incorporate sufficient fault tolerance properties so that they can 


accurately detect an envelope crossing and adjust the control settings accordingly. 
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APPENDIX 


COMPUTER PROGRAMS 


6THIS IS PROGRAM FORCER.M. IT COMPUTES THE BIASED 
BIFURCATION CUSP FOR 

S6DELTA SATURATED INCLUDING DRAG TERMS. IT USES THE FZERO 
FUNCTION TO 

S6SOLVE THE EQUATIONS. 


global W Bl Un Zw Zdlt A XA CD rho F xf; 
SESTABLISH GEOMETRIC PARAMETERS 


W: 5'2556.7363 
percent=input(’Enter percent of excess buoyancy '); 


Bl = (1+percent/100) *W; 

F = input(’Enter vertical force ‘'); 

xf = input(’Enter distance from amidships where force 
acts’); 


A = 19.8473; 
KA. = 0.20126; 
CD=input (‘Enter CD ‘}; 


Ly =. 561.327 
Gg = 32.2; 

m = W/g; 

tho = 1.94; 
eS 1S 192; 
<b] -0- 
Alphas = +1; 
Alphab = Q; 
2g = 0.1; 

ZB = 0% 


SNON-DIMENSIONING FACTORS 


Nd. = -0.5*rne*L" 2: 

naz ‘= 0.5*rne*l°3> 

na. =" °045"*rno* GL" 4+ 

na4. =. 0.5*7no*h"5; 
SHYDRODYNAMIC COEFFICIENTS 
Zzqdnd = -6.33e-4; 

Z2wand = -1.4529e-2; 

Zqnd = 7.545e-3; 

Zzwnd = -1.391e-2; 
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zZzas = -5.603e-3; 

Zep: = 0..5720s- 

Zdltnd = (Alphas*Zds+Alphab*Zdb) ; 
Mqdnd = -8.8e-4; 

Mwdnd= -5.61e-4; 

Mgnd = -3.702e-3; 

Mwnd = 1.0324e-2; 

Mds = -0.002409; 

Mdb = 0.5*(-Mds) ; 

Mditnd = (Alphas*Mds+Alphab*Mdb) ; 


zqd nas *Zqqnds 
Zwd nd2*Zwdnd; 
Zq = nd2*Zqgnd; 

aw = ndil*Zwnd; 
ZaQ1t = ndil*Zdltnd; 


Mgd = nd4*Mgdnd; 
= nd3*Mwdnd; 
Mg = nd3*Mgqnd; 
Mw = nd2*Mwnd; 
Malt = nd2* Mdlitnd; 


SESTABLISH SPEED RANGE 


Ga -0iy 3000522 ..6> 
Fn = U./sqrt(g.*zg) ; 
for n=l*length (Vv); 


SSOLVES FORCE EQUATION. USING FZERO FUNCTION 

th(n)= £zero(’gsatur’, (-Zdlt/Zw) ,0.00001) ; 

thi(n)= fzero(’gqsaturl’, (-Zdlt/Zw) ,0.00001); 

SCOMPUTES RELATIONSHIP BETWEEN U AND Xgb WITH RESULTS FROM 
ABOVE EON 


xg(n) = (Mw*U(n)%*2*tan(th(n))-zg*W*sin(th(n))- 
(05:5) =ENO"CD* A" XA*U (hn) 2 ww 
tan(th(n))*abs(tan(th(n)))+Md1lt*U(n)%*2* (0.4)- 
F*xt)/ (W*cos(th(n))).; 

end; 


XobL. = xg/L: 
Lp = sqrt (U.%*2/3.22) ; 


for n=1; length (Vy); 


xgl(n) = (Mw*U(n)*2*tan(thi1(n))-zg*W*sin(thi1 (n))- 
(0.5) *rno*CDT At KA « os 
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U(n)*2*tan(thi(n))*abs(tan(thi (n)))+Md1lt*U(n) *2* (-0.4)- 
BE 7 aces 

(W*cos(thi(n))); 
end; 


XobL1l= xgl1/L; 

plot (XgbL,Fn,XgbL1,Fn),grid 
title(’Saturation Cusp ’); 
xlabel (‘Xgb’); 

ylabel (’Fn’'); 


STHIS IS PROGRAM GSATUR.M. IT ESTABLISHES THE FUNCTION TO 
BE SOLVED IN 
S6THE BIFSU.M PROGRAM. 


function y = gsatur(th); 

global W Bl Un Zw Zdlt A XA CD rho F xf; 

y= Zw*U(n)*2*tan(th)+(W-B1)*cos(th) + Zdlt*U(n)*2* (0.4)... 
~(0.5) *rho*CD*A*U(n)*2*tan(th) *abs(tan(th)) - F; 


STHIS IS PROGRAM GSATUR1.M. IT ESTABLISHES THE FUNCTION TO 
BE SOLVED IN 
S6THE BIFSU.M PROGRAM. 


function. vl = osaturl(chl): 

global W Bl Un Zw Zdlt A XA CD rho F xf; 

yl =Zw*U (n) *2*tan(th1)+(W-B1)*cos(th1) + zZd1lt*U(n)*2* (- 
Ores aera 

-(0.5) *rho*CD*A*U(n)%*2*tan(thi1)*abs(tan(thi)) - F; 
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